The Mid-atlantic Ridge (31°S–34°30′S): Temporal and spatial variations of accretionary processes     

P. J. Fox  N. R. Grindlay  K. C. MacDonald 《Marine Geophysical Researches》1991,13(1):1-20

The ridge located between 31° S and 34°30′S is spreading at a rate of 35 mm yr⁻¹, a transitional velocity between the very slow (≤20 mm yr⁻¹) opening rates of the North Atlantic and Southwest Indian Oceans, and the intermediate rates (60 mm yr⁻¹) of the northern limb of the East Pacific Rise, and the Galapagos and Juan de Fuca Ridges. A synthesis of multi-narrow beam, magnetics and gravity data document that in this area the ridge represents a dynamically evolving system. Here the ridge is partitioned into an ensemble of six distinct segments of variable lengths (12 to 100 km) by two transform faults (first-order discontinuities) and three small offset (< 30 km) discontinuities (second-order discontinuities) that behave non-rigidly creating complex and heterogeneous morphotectonic patterns that are not parallel to flow lines. The offset magnitudes of both the first and second-order discontinuities change in response to differential asymmetric spreading. In addition, along the fossil trace of second-order discontinuities, the lengths of abyssal hills located to either side of a discordant zone are observed to lengthen and shorten creating a saw-toothed pattern. Although the spreading rate remains the same along the length of the ridge studied, the morphology of the spreading segments varies from a deep median valley with characteristics analogous to the rift segments of the North Atlantic to a gently rifted axial bulge that is indistinguishable from the shape and relief of the intermediate rate spreading centers of the East Pacific Rise (i.e., 21°N). Like other carefully surveyed ridge segments at slow and fast rates of accretion, the along-axis profiles of each ridge segment are distinctly convex upwards, and exhibit along-strike changes in relief of 500m to 1500 between the shallowest portion of the segment (approximate center) and the segment ends. Such spatial variations create marked along-axis changes in the morphology and relief of each segment. A relatively low mantle Bouguer anomaly is known to be associated with the ridge segment characterized by a gently rifted axial bulge and is interpreted to indicate the presence of focused mantle upwelling (Kuo and Forsyth, 1988). Moreover, the terrain at the ends of each segment are known to be highly magnetized compared to the centers of each segment (Carbotte et al, 1990). Taken together, these data clearly establish that these profound spatial variations in ridge segment properties between adjoining segments, and along and across each segment, indicate that the upper mantle processes responsible for the formation of this contrasting architecture are not solely related to passive upwelling of the asthenosphere beneath the ridge axis. Rather, there must be differences in the thermal and mechanical structure of the crust and upper mantle between and along the ridge segments to explain these spatial variations in axial topography, crustal structure and magnetization. These results are consistent with the results of investigations from other parts of the ridge and suggest that the emplacement of magma is highly focused along segments and positioned beneath the depth minimum of a given segment. The profound differences between segments indicate that the processes governing the behavior of upwelling mantle are decoupled and the variations in the patterns of axis flanking morphology and rate of accretion indicate that processes controlling upwelling and melt production vary markedly in time as well. At this spreading rate and in this area, the accretionary processes are clearly three-dimensional. In addition, the morphology of a ridge segment is not governed so much by opening rate as by the thermal structure of the mantle which underlies the segment.  相似文献   

Comparison of methods to model the gravitational gradients from topographic data bases   总被引：7，自引：0，他引：7  

Christopher Jekeli  Lizhi Zhu 《Geophysical Journal International》2006,166(3):999-1014

A number of methods have been developed over the last few decades to model the gravitational gradients using digital elevation data. All methods are based on second-order derivatives of the Newtonian mass integral for the gravitational potential. Foremost are algorithms that divide the topographic masses into prisms or more general polyhedra and sum the corresponding gradient contributions. Other methods are designed for computational speed and make use of the fast Fourier transform (FFT), require a regular rectangular grid of data, and yield gradients on the entire grid, but only at constant altitude. We add to these the ordinary numerical integration (in horizontal coordinates) of the gradient integrals. In total we compare two prism, two FFT and two ordinary numerical integration methods using 1" elevation data in two topographic regimes (rough and moderate terrain). Prism methods depend on the type of finite elements that are generated with the elevation data; in particular, alternative triangulations can yield significant differences in the gradients (up to tens of Eötvös). The FFT methods depend on a series development of the topographic heights, requiring terms up to 14th order in rough terrain; and, one popular method has significant bias errors (e.g. 13 Eötvös in the vertical–vertical gradient) embedded in its practical realization. The straightforward numerical integrations, whether on a rectangular or triangulated grid, yield sub-Eötvös differences in the gradients when compared to the other methods (except near the edges of the integration area) and they are as efficient computationally as the finite element methods.  相似文献   

Rescaled-range and power spectrum analyses on well-logging data   总被引：6，自引：0，他引：6  

Chun-Feng Li 《Geophysical Journal International》2003,153(1):201-212

  相似文献   

SH Wave Propagation in Viscoelastic Media     

Daley  P.F.  Krebes  E.S. 《Studia Geophysica et Geodaetica》2004,48(3):563-587

When comparing solutions for the propagation of SH waves in plane parallel layered elastic and viscoelastic (anelastic) media, one of the first things that becomes apparent is that in the elastic case the location of the saddle points required to obtain a high frequency approximation are located on the real p axis. This is true of the branch points also. In a viscoelastic medium this is not typical. The saddle point corresponding to an arrival lies in the first quadrant of the complex p-plane as do the branch points. Additionally, in the elastic case the saddle point and branch points lie on a straight line drawn through the origin (the positive real axis in the complex p-plane), while in the viscoelastic case this is generally not the case and the saddle point and branch points lie in such a manner as to indicate the degree of their complex values.In this paper simple SH reflected and transmitted particle displacement arrivals due to a point torque source at the surface in a viscoelastic medium composed of a layer over a half space will be considered. The path of steepest descent defining the saddle point in the first quadrant will be parameterized in terms of a real variable and the high frequency solutions and intermediate analytic results obtained will be used to formulate more specific constraints and observations regarding saddle point location relative to branch point locations in the complex p-plane.As saddle point determination for an arrival is, in general, the solution of a non-linear equation in two unknowns (the real and imaginary parts of the complex saddle point p ₀), which must be solved numerically, the use of analytical methods for investigating this problem type is somewhat limited.Numerical experimentation using well documented solution methods, such as Newton's method, was undertaken and some observations were made. Although fairly basic, they did provide for the design of algorithms for the computation of synthetic traces that displayed more efficient convergence and accuracy than those previously employed. This was the primary motivation for this work and the results from the SH problem may be used with minimal modifications to address the more complicated subject of coupled P-SV wave propagation in viscoelastic media.Another reason for revisiting a problem that has received some attention in the literature was to approach it in a fairly comprehensive manner so that a number of specific observations may be made regarding the location of the saddle point in the complex p-plane and to incorporate these into computer software. These have been found to result in more efficient algorithms for the SH wave propagation and a significant enhancement of the comparable software in the P-SV problem.  相似文献   

An analytical solution for describing the transient temperature distribution in an aquifer thermal energy storage system     

Kuang‐Yi Li  Shaw‐Yang Yang  Hund‐Der Yeh 《水文研究》2010,24(25):3676-3688

A mathematical model is developed for predicting the temperature distribution in an aquifer thermal energy storage (ATES) system, which consists of a confined aquifer bounded from above and below by the rocks of different geological properties. The main transfer processes of heat include the conduction and advection in the aquifer and the conduction in the rocks. The semi‐analytical solution in dimensionless form for the model is developed by Laplace transforms and its corresponding time‐domain solution is evaluated by the modified Crump method. Field geothermal property data are used to simulate the temperature distribution in an ATES system. The results show that the heat transfer in the aquifer is fast and has a vast effect on the vicinity of the wellbore. However, the aquifer temperature decreases with increasing radial and vertical distances. The temperature in the aquifer may be overestimated when ignoring the effect of thermal conductivity. The temperature distribution in an ATES system depends on the vertical thermal conduction in the rocks and the horizontal advection and thermal conduction in the aquifer. The present solution is useful in designing and simulating the heat injection facility in the ATES systems. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

Interpretation of aeromagnetic anomalies over part of Southeastern Nigeria using three-dimensional Hilbert transformation     

C. O. Ofoegbu  N. L. Mohan 《Pure and Applied Geophysics》1990,134(1):13-29

We present the results of the application of three-dimensional Hilbert transformation to the analysis of airborne total field magnetic anomalies over part of Southeastern Nigeria. This study not only substantiates the usefulness of 3-D Hilbert transforms in the interpretation of magnetic anomaly maps but also more clearly delineates the structural pattern of the area, of study. Results from the previous study are discussed in relation to the results of previous geological and geophysical studies of the area.  相似文献   

An examination of the use of the logratio transformation for the testing of endmember hypotheses     

Ross M. Renner 《Mathematical Geology》1991,23(4):549-563

Given a compositional dataset in the absence of any prior information on any mixing process which may have formed it, a complete analysis of mixtures determines three distinct types of estimates in order. These are: (i) the estimate of the number of endmembers or fixed source compositions, of which all the sample compositions of the dataset must be approximate mixtures; (ii) the estimated compositions for each of these chosen number of endmembers; and (iii) the estimated contributions of each of these endmember estimates to each sample. Traditionally, the estimate for the number of endmembers has been assessed either by mapping or by inspection of the coefficients of determination between the observed and estimated variables. Mapping entails the plotting on a map of the region from which the samples were taken, either the contours of the contributions of each endmember to each sample, or some other portrayal of the distribution of endmember abundances. Because it requires the complete analysis, assessment by this method is too elaborate except for final confirmation and display. Alternatively, choosing a number of endmembers, which result in suitability high coefficients of determination for all or most variables, may account for elements which are not part of the conjectured mixing process or, worse, may result in the identification of endmembers which may never in fact have existed. Such an error is similar to overspecifying a multiple regression model. So, the obvious starting point from which to assess the validity, or otherwise choice of endmember numbers, is to examine the matrix of residuals. The differences between the logratio-transformed observed and estimated data form an array of residual logratios. A linear combination of these may be formed for each sample, which, under a random perturbation assumption, should follow a univariate normal distribution. Whether or not this scalar is normal can be readily tested. It can also be examined graphically for such desirable qualities as symmetry when the test for normality may be too severe. This procedure is employed to assess the decompositions of the U.S.G.S. Mid-Pacific data and the Nazca Plate Surface sediments.This paper was presented at the 18th Geochautauqua, Newark, Delaware, 13–14 October 1989.  相似文献   

一类广义积分的算法及实现     

陶诏灵 《南京气象学院学报》2002,25(1):100-104

Laplace变换与B样条在社会生活的各个领域都有着广泛的应用，考虑将二者联系起来，运用B样条方法求解Laplace积分。对Laplace积分具有导数不连续点的像原函数，运用重节点B样条进行拟合，再引入差商，以两个不同出发点导出Laplace积分的算法。并在均匀网格分划下，给出一系列解析公式，然后用算例编程计算，实例表明简单、快速而且便于计算机实现。  相似文献   

半空间饱和土在内部简谐水平力作用下的Ge函数   总被引：6，自引：0，他引：6       下载免费PDF全文

陆建飞  王建华  沈为平 《地震学报》2001,23(2):185-191

根据积分变换方法得出了半空间内部作用简谐水平力时的Gree函数．首先，利用Hankel积分变换方法，直接对频域内的Biot波动方程进行求解，得出Biot波动方程的通解；利用通解和半空间内部作用水平力时边界上的边界条件，以及力作用面上的连续性条件，可以得出上述边值问题的解；对于边值问题在变换域内的解进行相应的逆变换，就可以得出频域内的Gree函数．本文得到的线弹性退化解与文献中的结果吻合．最后，文中给出了两个算例．   相似文献   

The Green''''s function of the harmonic horizontal force applied at the interior of the saturated half-space soil     

陆建飞  王建华  沈为平 《地震学报(英文版)》2001,(2)

IntroductionThe wave propagation problems in saturated soil are very important for the civil engineering, geophysics and seismology. Biot (1956,1962) established the theory of wave propagation in saturated soil firstly, and hereafter many researchers have used Biot theory to study wave propagation problems in saturated soil. By using integral transform and potential function method, Philippacopoulos (1988) studied the Lamb(s problem of a vertical point force applied to the surface of saturate…  相似文献